Correlations in the Adiabatic Response of Chaotic Systems 
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Adiabatic variation of the parameters of a chaotic system results in a fluctuating reaction force. 
In the leading order in the adiabaticity parameter, a dissipative force, that is present in classical 
mechanics was found to vanish in quantum mechanics. On the time scale t, this force is proportional 
to I{t), the integral of the force- force correlation function over time t. In order to understand the 
crossover between the classical and the quantum mechanical behavior we calculated I(t) in random 
matrix theory. We found that for systems belonging to the Gaussian unitary ensemble this crossover 
takes place at a characteristic time (proportional to the Heisenberg time) and for longer times I{t) 
practically vanishes, resulting in vanishing dissipation. For systems belonging to the Gaussian 
orthogonal ensemble I{t) ~ 1/t and there is no such characteristic time. I{t) is calculated for 
various models and the relation to experiment is discussed. 



PACS numbers: 05.45.Mt, 05.45.Ac, 05.45.-a 

Friction is usually described in statistical physics as 
transfer of energy to a system that consists of an infinite 
number of degrees of freedom in an irreversible way. The 
irreversibility is a result of the complexity of the motion 
of the many body system. A natural question is whether 
the coupling to a small chaotic system can exhibit dissi- 
pation resulting from its complex dynamics. Mechanisms 
of this nature were first introduced in the context of nu- 
clear physics llj , and various aspects of it were discussed 
in the past P-pl . In the present work we study, following 
Berry and Robbins (BR) Q and Jarzynski [^, a specific 
model where the chaotic system is defined by the Hamil- 
tonian H (R, z) where z = (p, r) are the phase space co- 
ordinates of the chaotic system and R is an adiabatically 
varying parameter. The system is chaotic for each value 
of R. The crucial feature of the system we studied in 
this work H is that it exhibits a wide separation of time 
scales - the evolution of the chaotic system, characterized 
by the time scale Tf^g-t, is so rapid that it explores all of 
the phase space accessible to it before the parameter R, 
characterized by the time scale Tgi^ , changes apprecia- 
bly. The adiabaticity parameter is e ~ ^fast/"^slow ^^^ 
average generalized force, that is applied by the chaotic 
system in question on the system that forces the adia- 
batic variation of R, is given by: 



r(ra)-- fdz p{z,Ta)dRniz,K{Ta)), 



(1) 



where p(z,Ta) is a normalized probability density in the 
fast particle phase space. If R is a position space co- 
ordinate, F(Ta) reduces to a regular force. It can be 
expanded in powers of the adiabaticity parameter e as 
was done by BR |B| and Jarzynski . The present work 
follows closely the formalism of BR. 

Berry and Robbins H were able to calculate the force 
acting on the slow particle up to first order in e: 



eFi 



(2) 



— 9fl;i?(R) where E (R) is chosen such that the phase 
space volume of the fast particle, il.{E{IV) ,R), is con- 
stant 1^. The leading correction to Fo includes a ve- 
locity dependent force Fi(Ta), and describes two quali- 
tatively different forces. The first of these is geometric 
magnetism, that is related to the Berry phase H. The 
second one is related to deterministic friction §,|l^,g . A 
central question is under which conditions friction due to 
the velocity dependent force Fi is found. 

For the exploration of the existence of friction it is 
sufficient to study the case where R is replaced by a 
scalar X. In this case the velocity dependent force Fi (ra) 
reduces to 

F,{Ta) = -X ^-^dE P I{t = ^)]e^e(X) . (3) 

where T.{E,X) = dE^{E,X) and: 



I{t) = / C{t')dt' 
Jo 



(4) 



The fiuctuating force-force correlation function at time 
difference t' is: 



where {■ ■ ■) ^ x denotes the microcanonical average over 
the energy surface with a fixed value of the parameter 
X. The volume enclosed by this surface is an adiabatic 
invariant ||^. The fiuctuation of the energy is 7Y (z, X) = 
7i(z,Ar) — E{X). A finite value of the integral /(oo) 
is required for friction. Within classical mechanics this 
integral is indeed positive. It was shown by BR that in 
quantum mechanics the dissipative part of the force Fi 
vanishes. To understand the reason for this discordance 
(g) was calculated and in quantum mechanics it takes the 
form: 



To leading order, the force is given by the classical 
analogue of the Born-Oppenheimer force FQ^{Ta) = 



C{t) 
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(6) 



where the Em are the eigenenergies of the chaotic sys- 
tem. The initial state is n and the dependence on it is 
not important for the present work. The sum is over 
TO. For this correlation function the integral /(oo) of (0) 
vanishes. In order to understand how the crossover be- 
tween the classical and the quantum behavior occurs, it 
is instructive to calculate the integral of the correlation 
function over a finite time. Taking the classical limit 
h ^' for any finite t and then the limit i — > oo should 
result in a non- vanishing value of /(oo), while for any 
finite value of h, /(oo) should vanish. The friction on 
the time scale t is proportional to I{t) as can easily be 
inferred from (^. For systems whose classical dynamics 
is chaotic, the energy levels are distributed according to 
random matrix theory (RMT) |lj]. The long time be- 
havior of I(t) is determined by the levels nearest to n, 
namely nil, as can be seen from (Jq). The natural ques- 
tion to ask is whether there is a characteristic time scale 
for the crossover between the quantum behavior of the 
integral I{t) and its classical behavior. The most naive 
answer to this question is that the characteristic time 
scale is the Heisenberg time. On the other hand, one can 
argue that there is no time scale for this crossover at all 
[ jl2[ . In RMT the probability for two consecutive levels 
to be separated by an energy difference s behaves like s^ 
for small spacings [|l3|,0. Consequently, {Ip{t)) ~ t~l^ 
for long times. Here (. . .) denotes the RMT ensemble 
average, and Ip is the integral (0) for some (3. The (. . .) 
will be dropped from Cp and Ip for notational simplicity. 

For the Gaussian orthogonal ensemble (GOE) (/3 = 1) 
one indeed finds that Ip{t) decays like 1/i, but for the 
Gaussian unitary ensemble (GUE) (/3 = 2) one finds that 
it decays like a Gaussian with a characteristic time pro- 
portional to the Heisenberg time or vanishes after the 
Heisenberg time depending on the parametric depen- 
dence on X. Why is the nature of the decay of Ipit) 
important? There is the quantum-classical discordance 
that has already been mentioned, and one would like to 
analyze the scale that is required to observe the crossover 
between the regimes. In addition, the model discussed in 
the present work is relevant for some experimental sit- 
uations. Consider for example a molecular beam pre- 
pared in a classical configuration, where initially many 
levels are substantially populated. The beam travels in 
a slowly varying field [l5[ . Consequently the internal dy- 
namics in the molecules is in a slowly varying potential. 
Another example is of quantum dots where parameters 
are varied adiabatically like in pumping experiments, but 
with closed dots, so that their spectrum is discrete [16|. 
A dramatic change in the energy absorption is predicted 
when a magnetic field is applied and the transition from 
GOE to GUE takes place. 

The time over which the correlation function decays 
should be compared with other time scales present in the 
specific system studied. One such time scale is T2 ~ e~^, 
which is the time scale for the breakdown of the first order 



of the multiple scale expansion in e. Non-perturbative ef- 
fects, such as Landau-Zener tunneling, become important 
on a time scale of T^z- In realistic experiments there is 
also the time scale for quantum decoherence T^. In order 
to observe the classical to quantum crossover discussed 
in the present work (Ipit)) should exhibit substantial de- 
cay for t ^ inm{T2,TLz) and of the order of T^. The 
energy absorption by the internal degrees of freedom is 
proportional to I(3{T^). 




FIG. 1. The integral of the correlation function for GOE. 
Numerical results for AT = 3 (o), A^ = 13 (v), N = 53 (d) 
and A*' = 103 (a) are shown. Also shown is the large N 
approximation (Eq. h4l) (line). 




FIG. 2. Similar to Fig. |l| for GUE, compared to the large 
N approximation (Eq. um. 

We shall use the Hamiltonian introduced by Austin 
and Wilkinson Q and model a parameter dependent sys- 
tem by the N x N random matrix: 

H{X)=HiCosX + H2SmX, (7) 

where Hi 2 are N x N random matrices from the same 



RMT ensemble. There are three advantages to working 
with H{X): (a) it belongs to the same ensemble that 
Hi 2 belong to; (b) the derivatives of its matrix elements 
belong to the same ensemble; (c) the matrices H{X) and 
H'{X) = dH{X)/dX are statistically independent. If we 
insert H{X) and H'{X) into Eq. ^ and then perform the 
ensemble average, we obtain: 




Cp{t) 



where H'{X)n,; 



H'{X), 



cos 



{En — E„ 



and 



(8) 



de- 



dH{X)/dX 

notes RMT ensemble averaging. 

The statistical independence of H'{X) and H{X) im 
plies 



Cp{t) = |3^i■ 



2 Y^ 




t 



{En — En 



(9) 



The fact that H'{X) belongs to the same ensemble as 
H{X), leading to {\H' {X)n,„^\' ) = ( li?™,™!') = /3m' 
for m ^ n^ was used. 

We would like to make the connection between 
Ci3{t)/ fifj? and the form factor: 



K{t) 



1 



r{E) 



(p{E + e/2p)p{E-e/2p))-l 



i27reT 



de, 



(10) 



where p{E) = J^i ^{Ei — E) is the density of states and 
'p{E) is the smoothed density of states. The variable e is 
the energy measured in units of the mean level spacing 
1/J){E) and r = t/Tn is time in units of the Heisenberg 
time, Th = h'p{E). In what follows units where Pp^Tn = 
1 will be used. Eq. M can be written in the following form: 



Cp{T) 



-p\E) 



p{E + e/2p)p{E-e/2p))~S{e) 



i2-KeT 



de. 



(11) 



Comparing the last equation with (10) one can see that: 

Cp{t)IPp^ = K{t) + S{t)-1. (12) 

In this work we arc mainly interested in the time integral 
of the correlation function (^) : 



Ip{T) = 1/2- 



I dT'(l-K{T'] 



(13) 



In the limit t ^ oo the RHS is just R2{e ~ 0), the two 
point spectral correlation function at zero energy separa- 
tion. It vanishes as a result level repulsion. 

In order to perform actual calculations we make use of 
the well known form factor for GOE and GUE [Q. It is 
standard to define 6(r) = 1 - K{t). For GOE (see 



p. 137): 6(r) = 1 - 2r + Tln6+ for r < 1 and 6(t) = 
— l-|-Tln[6+/6_] for T > 1 where b± = 2t±1, leading to: 




i(r2-i)ln[l + 2r]] t<1 



l{l-r) + \{r--\)\n 



For r — > oo it falls off asymptotically as 
/i(t) ^ l/12r. 



2t + 1 



T> 1 



(14) 



(15) 



For GUE (see |lj p. 95): 6(r) = 1 - t for r < 1 and it 
vanishes for r > 1, from which one obtains: 



h{r) = 



_ ri/2-[r-rV2] 







T < 1 
T> 1 



(16) 



The behavior for the Gaussian symplectic ensemble 
(GSE) is very similar. 

In order to compare the analytical results that hold in 
the infinite N limit with numerical data, ensembles of 
N X N matrices Hi and H2 of (|^), belonging to GOE 
or GUE were generated numerically and the results are 
presented in Figs. & §. 

The model (M) is very specific in its dependence on the 
parameter X. An important property of this model is 
the statistical independence between H{X) and H'{X). 
Such independence holds to a good approximation for 
disordered systems [n7|. It is reasonable to make this 
approximation also for RMT models of chaotic systems. 
The reason is that most eigenstates look random, are 
statistically independent of the eigenvalues and therefore 
for many types of perturbations the matrix elements of 
H'{X) will look random and independent of the spec- 
trum. Although this argument is reasonable for many 
types of parametric dependence it is clearly not general. 
For the asymptotic behavior much less is required, since 
the long time asymptotics is dominated by the nearest 
neighbor level spacing. The reason for this dominance is 
that if r ^ 1 the terms in the sum (||) oscillate wildly as 
a function of m, so that the important net contribution 
is from the terms nearest to being stationary. These are 
obviously m = 71 ± 1. Introducing also the crucial as- 



sumption that that fiuctuations of 
be ignored, we find that in RMT: 



{H'{XK,n 



±1 



can 



Cp{r)/Pp^ 



/•CXD 

2 I ds Pp{s) cos {27TTS), (17) 

Jo 



where s is the nearest neighbor level spacing in units of 
the mean level spacing AE — l/p(0) and Pf3{s) is the 
distribution of s. The integral of the correlation function 
is: 



Ipir) 



1 rds^^sm{2TTTs). 

T^ Jo S 



(18) 



For the nearest neighbor level spacing distribution we use 
the Wigner surmise (see Eq. 202 in |l3]): 



Puis) 



c s' 



' exp [- 



as 



(19) 



with a = 7r/4 and c = 7r/2 for GOE, and a = A/tt and c = 
32/7r2 for GUE. The integral ^) can be calculated for 
these distributions. For GOE one finds that for large t it 
takes the form: /i(t) ~ 1/Attt, that is extremely close to 
(|l|). For GUE one finds: ^2(t) = 2r exp [-tt^t^/A] , that 
for large r is very close to (|6|). Here too the behavior 
for GSE is very similar to that of GUE. 

We found that for the RMT ensembles /^(r) of (|) is 
dominated by the nearest neighbor level spacings result- 
ing in the approximation dlH) for //3(t). Under the as- 
sumptions leading to this approximation I /sir) can be cal- 
culated for the nearest neighbor level spacing distribution 
( p^ ) with arbitrary /3, a and c even if H does not belong 
to an invariant Gaussian RMT ensemble. The integral 
of the correlation function (llg) takes the form //^(t) — 



0-1 



exp [—s\ smsy 



(c/7ra^/2) Ipiy) where X^(y) = J^s s 
with y = 2'KT/^/a. It satisfies the ordinary ordinary dif- 
ferential equation: 



y 



P. 



^^/3(y) + n-Mv) + ^My) = 



2dy 



with the boundary conditions I^ (0) — and dIi3{Q)/dy — 
T[{j3 + I) /2]/2. Making the substitution l^iy) = 
//3(j//-\/2) exp [— 2/^/8] and changing to the variable x = 
?//\/2, one arrives at a new differential equation for / that 
is a well known equation, and its solutions are Parabolic 
Cylinder Functions @: U [1/2 - /3, x] , V [1/2 - /?, x]. 
For arbitrary (3, its solution is a linear combination of 
these two functions. From the initial conditions and the 
known behavior of U and y at and oo the asymptotic 
behavior of Ipir) is found to be of two types. For (3 ^ 2n 
(n = 1,2,3,...): //3(t) ^ r^^, as r ^ oo, while for 
j3^2n{n^ 1,2,3,...): /^(r) - exp [-TrVVajr^-^, as 
T — !■ cx). This is precisely the type of behavior found for 
the RMT ensembles treated explicitly. 

It was shown that for the model ^ there is a big qual- 
itative difference between the integrals of the correlation 
functions of the RMT ensembles (Eqs. |l| & ^. It was 
argued that the qualitative difference holds for a wide va- 
riety of models. For models where the mean level spacing 
is ([l9|), //3(r) decays like the power-law t"'', except for (3 
that are positive even integers, for which it decays like a 
Gaussian with a characteristic time that is proportional 
to the Heisenberg time. Only for even values of (3 the in- 
tegral (ttq) can be written as an integral over an analytic 
function in the interval [— oo, oo], and is dominated by a 
saddle point in the complex plane. 
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